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INTRODUCTION 


There  are  at  least  three  examples  of  multiple-crack  patterns  which  exhibit  some  degree 
of  stability.  The  first,  which  exhibits  an  extremely  high  degree  of  stability,  concerns 
collinear  periodic  arrays  of  cracks  emanating  from  a  pre-existing  array  of  collinear, 
equally  spaced  rivet  holes.  Because  of  the  principal  driving  force  (the  force  exerted  by 
the  rivet  itself)  and  the  known,  analytic,  value  of  stress  intensity  for  such  configurations 
there  is  an  inexorable,  self-correcting,  drive  to  create  multiple,  equaHength  cracks  In 
such  a  configuration.  Since  the  force  exerted  by  the  rivet  is  automatically  adjusted  (via 
a  change  in  local  compliance)  and  because  the  crack  tip  nearest  any  rivet  has  a  higher 
stress  intensity  than  the  more  distant  tip  [1]  a  symmetric,  equaNength  array  is  ensured. 
The  second  example,  which  exhibits  somewhat  less  stability,  concerns  arrays  of  radial, 
periodic  cracks  emanating  from  a  circular  boundary  with  an  applied  displacement 
boundary  condition,  for  example  cracks  in  glass  after  the  passage  of  a  bullet,  cracks  in 
ice  following  a  pointed  indentation  and  cracks  in  dried-out  tree  stumps. 

A  third  example,  which  exhibits  moderate  stability  and  is  considered  herein,  concerns 
multiple-edge  cracks  emanating  from  a  free  surface.  Multiple  cracking  from  the  free 
surface  of  an  (effectively)  semi-infinite  body  due  to  tensile  stresses  in  the  layer  at  (or 
near)  the  free  surface  is  a  natural  phenomenon.  The  surface  topography  of  the 
cracking  of  ice-wedge  polygons  in  Arctic  permafrost  [2],  of  mud  flats  in  Death  Valley, 
Figure  1(a),  and  of  craze-cracks  (heat-checking)  at  the  bore  of  a  gun  tube.  Figure  1(b), 
[3]  are  all  strikingly  similar,  yet  they  span  five  orders  of  magnitude  in  scale,  with  the 
maximum  plate  dimensions  for  ice  and  mud  being,  respectively,  22  metres  and  0.25 
metres  and  with  the  minimum  plate  size  for  craze-cracking  being  0.2mm  [3], 

Consider  a  system  of  multiple,  equal-length  edge  cracks  of  depth  a  and  spacing  2h. 

This  geometry,  together  with  some  relevant  stress  intensity  factor  values  [4],  is 
presented  in  Figure  2.  Lachenbruch,  in  his  lengthy  treatise  [2],  has  shown  that  linear 
elastic  fracture  mechanics  can  be  applied  to  edge-cracking  in  ice;  he  also  presented 
fracture  toughness  values  for  such  cracking.  However  Lachenbruch  did  not  explicitly 
consider  an  array  of  edge  cracks;  instead  he  calculated  contours  of  elastic  stress-relief 
at  the  surface  resulting  from  the  introduction  of  a  single  crack  of  depth  a  into  a 
pre-existing  tensile  stress  system  consisting  of  constant  tensile  stress  over  some  or  all 
of  the  prospective  crack  surface,  with  zero  stress  over  the  remainder.  He  concludes  that 
the  surface  stresses  resume  95%  of  their  pre-crack  level  at  a  distance  from  the  crack 
equal  to  9  x  crack  spacing  (i.e.  h/a  =  4.5).  Crucially  Lachenbruch  also  notes  that  a 
typical  observed  h/a  ratio  for  permafrost  cracking  is  1 .5. 

Bazant  and  co-workers  have  produced  several  papers  [5  -9]  in  which  are  analyzed  the 
sudden  cooling  of  a  rock  surface  and  the  associated  cracking  which  this  produces. 
Bazant's  analysis  is  based  upon  an  LEFM  study  of  the  energy  of  individual  cracks  and 
of  all  cracks.  Bazant  argues  that  crack-shedding  will  occur  as  a  proportion  of  the 
numerous  initial  equal-length  cracks  move  ahead  faster  in  order  that  the  system  may 
maximise  energy  release.  Crucially  Bazant  assumes  that  every  other  crack  is  shed 
(Shedding  ratio,  S  =  1/2)  and  calculates  the  crack  depth  at  which  shedding  occurs  as  a 
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function  of  the  driving  (residual  tensile  stress)  profile  and  concludes  that  shedding 
begins  at  h/a  ratios  which  depend  upon  the  local  residual  tensile  stress  and  are 
generally  in  the  range  0.166  <  hla  <  0.666.  Unfortunately  no  experimental  evidence  is 
presented  in  support  of  these  predictions. 

Both  Lachenbruch  (ice)  and  Bazant  (rock)  assume  that  the  progression  of  the  crack  is 
via  a  series  of  fast  fracture  events  when  crack  tip  stress  intensity,  K,  reaches  the 
Fracture  Toughness,  K^,  of  the  material.  However  in  the  case  of  a  gun  tube  it  is  known 
that  the  history  of  crack  growth  may  encompass  some  or  all  of:  Fast  Fracture; 
Environmental  Crack  Growth;  Fatigue  Crack  Growth.  Furthermore  there  is  no  longer  a 
single  driving  mechanism  since  global  residual  stresses  may  be  present  in  the  tube 
from  the  outset  (as  a  result  of  the  autofrettage  process)  and  others  are  induced  by 
thermal  and/or  phase  change  effects  near  the  bore  during  firing  of  the  first  few  rounds. 
Finally  the  cyclic  pressurization  of  the  tube  will  subsequently  serve  to  drive  cracks  by 
fatigue  crack  growth. 

Detailed  work  relating  to  the  gun  tube  problem  will  be  reported  elsewhere  [10].  However 
it  is  important  to  emphasize  that  all  of  the  above  modes  of  multiple-crack  extension  are 
K  driven,  so  that  any  K-based  crack-shedding  model  is  likely  to  be  applicable  across 
modes  of  crack  extension,  and  in  the  transition  between  modes,  when  these  occur 
consecutively  within  a  structure  as  in  the  case  of  a  gun  tube. 

Figure  3  illustrates  the  development  of  crack-shedding  in  the  case  where  two  out  of 
every  three  cracks  are  shed,  i.e.  the  shedding  factor  S  is  2/3.  Obviously  the  stress 
intensity  factor  and  associated  energy  changes  in  this  case  will  differ  from  Bazant's 
assumption  of  S=1/2.  It  may  therefore  be  useful  to  take  a  simplistic  approach  which 
permits  the  study  of  S  over  a  wide  range.  Re-plotting  the  results  presented  in  Figure  2 
on  a  logarithmic  scale  as  a  direct  function  of  h/a  we  obtain  the  K  variation  represented 
in  Figure  4  by  triangles. 

Note  that  the  variation  in  K,  although  normalized  by  cjh  does  properly  represent  the 
relative  variation  of  real  K  provided  h  is  held  constant  and  a  is  permitted  to  vary.  Figure 
4  also  shows  the  proportionate  variation  of  per  tip  energy  release  rate,  obtained  by 
squaring  K.  Finally  the  proportionate  variation  of  total  energy  release  rate  for  all  tips, 
obtained  by  summing  over  all  cracks  for  a  given  h/a,  is  also  plotted.  From  this 
presentation  we  conclude: 

a.  If  maximum  total  energy  release  rate  drives  crack-shedding,  either  this  is 
maximized  by  reducing  h/a,  an  effect  which  is  not  observed,  or  cracks  do  not 
shed,  which  also  contravenes  experimental  observation.  Total  energy  release 
rate  is  therefore  rejected  as  an  explanation  for  crack-shedding. 

b.  That  the  energy  release  rate  for  an  individual  crack  is  maximized,  with  fixed  h, 
by  increasing  h/a  to  1.5;  any  further  increase  in  a  will  actually  reduce  energy 
release  rate  slightly. 


2 


Having  rejected  total  energy  now  consider  the  option  of  crack-shedding  in  order  to 
maximize  the  energy  release  rate  of  an  individual  crack;  since  energy  release  is  directly 
proportional  to  K^,  all  subsequent  argument  is  based  directly  upon  K  variation.  Figure  5 
shows  proportionate  variations  in  real  K  values;  this  is  because  the  normalizing  term, 

Ko  =  c/hi ,  where  hi  is  the  initial  crack  spacing.  In  the  case  of  fixed  h  and  varying  a, 
movement  from  right  to  left  along  the  lower(continuous)  curve  represents  crack 
deepening,  which  occurs  incrementally. 

Starting  with  h  fixed,  from  a  typical  observed  initial  value  of  h/a  >  4,  K  increases  with 
increasing  a  until  at  h/a  =  1.5  (i.e.  2h/a  =  3).  K  cannot  increase  further  by  increasing  a 
and  the  only  option  is  to  begin  to  'shed'  a  proportion  of  the  cracks.  If  it  were  possible  to 
achieve  such  shedding  by  a  'healing'  process  in  which  a  proportion  of  the  cracks  have 
their  surfaces  brought  back  into  contact,  the  change  in  K  would  be  that  illustrated  for 
fixed  a  and  varying  h  in  the  upper  (dotted)  curve.  The  h/a  =  1 .5  point  on  the  lower  curve 
has  its  h/a  value  increased  by  a  factor  of  1/(1  -  S)  in  moving  to  hi a  =  1 .5/(1  -  S)  if  a 
proportion  S  of  the  cracks  were  'healed'  with  its  associated  stress  intensity  from  the 
lower  curve  increased  by  a  factor  of  7 1/(1  -S) .  'Healing'  is  clearly  physically  unreal  but 
will  serve  as  an  indicator  and  a  bound  for  the  shedding  process.  In  such  a  process  we 
note  that  K  would  increase  dramatically.  Since  only  integer  numbers  of  cracks  can  be 
shed,  and  periodic  symmetry  must  be  maintained,  a  crack  starting  from  h/a  =1.5  could 

move  to  h/a  =  3,  4.5,  6,  7.5 . etc..  These  integer  options  are  indicated  by  open  boxes 

along  the  dotted  curve.  Of  these  options  h/a  =  3  will  serve  to  achieve  a  38%  increase  in 
K,  4.5  gives  46%,  6.0  gives  47%  etc..  However,  the  h/a  =  4.5  option  also  permits  a 
transition  to  significant  amounts  of  energy  release  during  subsequent  crack  deepening 
and  is  a  prime  candidate. 

Returning  to  the  physical  unreality  of  'healing',  in  order  to  identify  a  more  physically 
acceptable  bound,  assume  that  in  the  natural  progression,  as  a  result  of  variations  in 
material  properties,  grain  size  etc.,  a  proportion,  S,  of  the  cracks  do  indeed  move  more 
slowly  than  their  neighbours  such  that  the  deeper  cracks  are  1.5  times  the  depth  of  the 
shorter.  A  heuristic  appeal  to  St  Venant  and  a  knowledge  of  the  behaviour  of  arrays  of 
deep  cracks,  indicates  that  the  influence  of  the  shorter  cracks  will  now  be  much  less 
significant.  (This  matter  is  discussed  in  more  detail  in  an  upcoming  section).  The  K 
variation  resulting  from  this  situation  is  shown  as  the  upper  discrete  points  (filled  boxes) 
in  Figure  5;  any  integer  point  on  the  'healing'  curve  has  its  h/a  ratio  decreased  by  1/3 
and  its  K  value  increased  in  proportion  to  the  ratio  of  the  K  values  relating  to  these  two 
h/a  values  extracted  from  the  lower  curve  in  moving  to  these  discrete  points.  The 
overall  argument  relating  to  the  relative  changes  in  K\s,  however,  unchanged. 

The  clear  implication  of  the  form  of  presentation  in  Figure  5  is  that,  to  increase  K, 
cracks  may  all  deepen  until  /7/a  =  1.5  followed  by  the  shedding  of  a  proportion  equal  to 
2/3  by  deepening  of  the  remainder  to  create  a  new  effective  h/a  =  4.5  whereupon  the 
option  of  increasing  fCby  further  shedding  is  removed  and  the  cycle  repeats  itself  by 
deepening  of  the  currently  active  (deeper)  cracks.  The  range  of  h/a  over  which 
deepening  and  shedding  progresses  is  therefore  strikingly  limited  to  1 .5  <  hla  <  4.5 
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which  is  precisely  the  region  within  which  the  short  and  deep  crack  limiting  solutions  do 
not  apply. 


EXAMPLE  AT  SMALL  SCALE  :  HEAT  CHECKING  OF  STEEL 

A  more  intuitive  presentation  which  permits  some  understanding  of  the  effect  of  stress 
gradients,  is  shown  in  Figure  6.  Typical  stresses  and  crack  depths  associated  with  the 
craze-cracking  of  the  surface  of  a  steel  gun  tube  are  assumed.  Two  prospective  stress 
profiles  are  shown  in  Figure  6(a).  They  consist  of  a  constant  pre-existing  stress  field  of 
SOOMPa  into  which  cracks  will  be  introduced  and  of  a  linearly  decreasing  field  with  a 
maximum  value  of  SOOMPa.  Stress  fields  which  increase  with  depth  are  not  observed  in 
practice. 

Accepted  approximations  to  stress  intensity  are  based  upon  direct  stress  normal  to  the 
prospective  crack  line  at  the  prospective  crack  tip  location,  designated  a*  [1 1, 12].  For  a 
single-edge  crack  the  variation  of  K  with  depth,  based  upon  AT  =  1.12a*  Jita  ,  [11],  is 
shown  as  a  heavy  curved  line;  Figure  6(b)  relates  to  constant  stress  and  6(c)  to  linearly 
reducing  stress. 

However  these  profiles  are  modified  in  the  case  of  multiple  cracking,  the  maximum 
value  at  small  h/a  being  effectively  'capped'  at  K  =  c  *  Jh  [^2].  Figures  6(b)  and  6(c) 
illustrate  such  capping  for  three  typical  h  values  {h  =  0.2mm,  0.6mm  and  1.8mm 
respectively). 

Consider  first  the  case  of  constant  stress.  Figure  6(b).  An  array  of  very  short  edge 
cracks  with  h  =  0.2mm  may  deepen,  by  increasing  a,  until  reaching  a  depth  of 
approximately  0.13mm  at  which  point  K  cannot  increase  further  without  the  shedding  of 
a  proportion  of  cracks  in  order  to  move  asymptotically  towards  the  heavy  (single-edge 
crack)  line.  Furthermore  such  a  shift  cannot  take  place  by  moving  vertically  from  a  point 
on  the  'capped'  horizontal  line,  since  this  is  equivalent  to  the  'healing'  of  a  proportion  of 
the  cracks.  A  schematic  of  such  a  shift  is  illustrated  as  is  the  case  of  subsequent 
capping  at  0.4mm  and  further  shedding.  In  the  case  of  constant  stress  these  shedding 
processes  may  be  expected  to  continue  to  repeat. 

Any  shift  may  be  expected  to  be  gradual  and  towards  a  point  which  is  outside  the 
'interaction  range',  a  term  coined  in  [13]  to  indicate  the  difference  in  depth  between  a 
longer  and  shorter  crack  within  a  periodic  array  of  such  cracks  at  which  the  influence  of 
the  shorter  crack  upon  the  stress  intensity  of  the  longer  crack  becomes  insignificant. 

The  determination  of  interaction  range  for  cracks  of  the  periodicity  and  within  the  h/a 
range  of  interest  here  is  not  addressed  in  [13],  A  discussion  relating  to  the  asymptotic 
bounds  on  'interaction  range’  is  in  preparation  [14].  In  order  to  proceed  with  the  present 
analysis  it  is  assumed  pragmatically  that  'interaction  range'  is  approximately  equal  to  h  . 
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In  the  nnore  physically  realistic  case  of  linearly  decreasing  stress,  Figure  6(c),  the  heavy 
curved  (single-edge  crack)  line  represents  an  upper  bound  on  possible  /<  values;  this 
curve  and  all  'capping'  lines  go  through  a  single  point  at  a  =  0.67mm  (the  point  of  zero 
stress).  The  argument  relating  to  the  point  at  which  shedding  begins  is  similar  to  the 
case  of  constant  stress,  but  the  availability  of  a  full  'interaction  range'  in  which 
shedding  can  occur  is  now  far  more  limited.  There  is  the  possibility  of  one,  or  at  most 
two,  shedding  cycles  before  crack  arrest  occurs.  Of  crucial  importance  is  the  threshold 
value  of  Kat  which  crack  extension  ceases.  This  may  relate  to  a  threshold  for  fast 
fracture,  fatigue  crack  growth  or  environmental  crack  growth.  The  latter  is  of  prime 
concern  in  the  early  stages  of  crack  formation  in  a  gun  tube  and  a  typical  value  of 
UMPaJin  (typical  of  hydrogen  cracking  in  gun  steel)  is  assumed.  With  this  additional 
bound  on  crack  extension,  designated  KSCC,  a  maximum  of  one  shedding  cycle  may 
occur  with  arrest  at  approximately  0.2mm  crack  depth. 

Note  that  the  use  of  o*  provides  an  acceptable  approximation  to  K  provided  the  stress 
gradient  is  not  excessive.  This  point  is  illustrated  in  Figure  7  which  repeats  the 
approximate  bounds,  based  upon  a*  from  Figure  6(c)  and  also  shows  calculations  of 
actual  K  profiles.  The  method  used  to  obtain  K  is  of  extremely  high  accuracy  (errors  < 
0.5%)  and  involves  the  modified  mapping  collocation  technique,  [15]  packaged  as 
weight  function  data  [16].  With  this  correction  to  K  profiles  the  predicted  crack  depth  at 
arrest  is  increased  slightly  to  approximately  0.22mm.  In  a  further  paper  [10]  a  more 
detailed  analysis  of  the  gun  tube  problem  is  undertaken  and  more  accurate  numerical 
methods  are  indeed  found  to  be  necessary  in  order  to  adequately  characterise  K  in 
such  high  stress  gradients. 


EXAMPLE  AT  LARGE  SCALE  ;  CRACKING  OF  ARCTIC  PERMAFROST 

To  illustrate  the  application  of  this  approach  at  large  scale  consider  the  problem  of 
cracking  of  Arctic  permafrost,  discussed  at  length  by  Lachenbruch  [2];  such  cracking  is 
reported  to  be  most  dramatic  in  early  February  each  year.  The  heavy  curve  in  Figure  8 
shows  the  variation  of  stress  as  a  function  of  depth  on  1  February  1954.  This  curve  is 
based  upon  temperature  variation  with  depth  taking  account  of  the  self-weight  of  the 
permafrost  and  is  extracted  directly  from  Figure  9  in  [2].  Figure  8  also  shows  the  value 
of  stress  intensity  for  a  single-edge  crack  and  for  arrays  of  cracks  based  upon  c*  and 
the  accepted  approximations  used  in  the  previous  section.  Three  examples  of  semi¬ 
crack  spacing,  h,  are  shown,  namely  h  =  1 ,  3  and  9  m. 

Interpreting  Figure  8  in  analogous  fashion  to  the  previous  section  multiple  cracks  with 
a<0.5m  and  with  h=1m  show  rapidly  increasing  K  until  'capped'  by  the  array  effect  at 
around  0.5m  depth.  K  may  then  only  increase  further  by  'shedding'  and  moving 
asymptotically  to  the  h=3m  curve. 

From  the  h=3m  curve  cracks  may  shed  to  h=9m  since  there  is  available  a  full 
'interaction  range'  (i.e.  3m)  for  shedding.  However,  once  established  at  h=9m  a  full 
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'interaction  range'  (now  9m)  is  not  available  for  shedding  and  the  cracks  are  effectively 
'trapped'  at  their  current  h  value  of  9m  and  will  be  forced  to  arrest  when  K  drops  to 
(approximately  0.3  MPam’'''  for  ice)  giving  a  crack  depth  of  around  5m  and  crack 
spacing  (2h),  or  plate  size,  of  18m.  These  figures  are  each  in  reasonable  agreement 
with  the  average  of  observations  reported  in  [2]. 

In  a  similar  fashion  to  the  preceding  section  more  accurate  (MMC)  K  calculations  were 
performed.  MMC  results  analogous  to  the  approximate  results  in  Figure  8  are  shown  in 
Figure  9.  The  way  in  which  crack-shedding  and  extension  occurs  is  clearly  similar  to 
that  argued  for  the  approximate  case,  but  with  crack  depths  extending  somewhat 
deeper  to  6.5m.This  figure  is  in  excellent  agreement  with  the  average  of  observations 
reported  in  [2]. 


EXPERIMENTAL  EVIDENCE 

Figure  10  shows  the  edge-cracking  which  developed  at  the  bore  of  a  large  caliber  gun 
tube  [3]  and  progressed  via  heat-checking  and  fatigue  processes,  both  X-driven.  The 
deepest  of  the  cracks  is  approximately  0.5mm  in  a  wall  of  50mm  total  thickness.  It 
appears  that  the  initial  cracking  exhibited  an  h/a  ratio  greater  than  4.  The  inferred  crack 
history  appears  to  follow  that  illustrated  schematically  in  Figure  3.  Crack  growth  clearly 
begins  at  h/a  of  approximately  4.5  and  continues  until  h/a  is  approximately  1.5, 
whereupon  shedding  occurs  and  the  process  is  repeated  twice  more.  The  ratio  S 
appears  to  be  2/3  throughout. 

If  such  shedding  does  occur  then  a  histogram  of  crack  depth  distributions  should 
appear  bimodal.  Gough  and  Morrison  [3]  have  counted  numbers  and  depths  of  cracks, 
and  at  least  two  of  the  associated  histograms  do  appear  bimodal. 

In  analyzing  the  surface  topography  of  cracked  concrete  Bazant  and  Cedolin  [9]  note 
that  an  energy-based  analysis  predicts  that  cracking  will  progress  by  the  associated 
surface  effect  of  the  creation  and  shedding  of  hexagonal  patterns;  at  each  shedding 
event,  hexagons  whose  'diameter'  is  three  times  larger  are  created.  Note  that  such  a 
surface  effect  is  equivalent  to  the  value  of  h  in  our  two-dimensional  analysis  increasing 
threefold  (i.e.  S=2/3),  as  predicted  from  the  simple  model  presented  herein.  Reference 
[9]  also  reports  that  such  hexagonal  cracking  is,  indeed,  observed  in  concrete.  It  is 
therefore  somewhat  surprising  that  throughout  References  [5  -  9]  a  value  of  S=1/2  is 
assumed. 

On  a  much  larger  scale,  Lachenbruch  [2]  reports  permafrost  ice  crack  h/a  ratios  of  1 .5 
wherein  h  may  be  as  large  as  22  metres. 

Finally,  recent  experimental  observations  on  ’island-delamination'  in  brittle  coatings  of 
varying  depth  [17]  produce  convincing  evidence  of  a  consistent  2h/a  ratio  of  3:1  for 
ceramic  films  of  up  to  0.4  mm  depth  and  paint  films  of  up  to  0.1  mm  depth.  In  the  terms 
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of  [17]  the  distance  2h  used  herein  is  equivalent  to  the  island  diameter  and  a  is 
equivalent  to  the  film  thickness  (i.e.  maximum  depth  of  cracks). 


CONCLUSIONS 

The  creation  and  subsequent  shedding  of  arrays  of  edge  cracks  is  a  natural 
phenomenon  ranging  from  heat-checking  of  gun  tubes  to  cracking  of  permafrost 
covering  five  orders  of  magnitude  in  crack  spacing.  The  shedding  behaviour  appears  to 
be  governed  by  energy  release  from  individual  cracks  rather  than  global  energy 
changes. 

The  model  developed  herein  predicts  that  all  cracks  will  deepen  until  a 
crack-spacing/crack-depth  ratio  {2h/a)  of  3.0  is  achieved,  at  which  stage  crack-shedding 
will  commence:  two  out  of  every  three  cracks  will  be  shed,  leading  to  a  larger  crack¬ 
spacing/crack-depth  ratio  at  which  stage  growth  of  all  currently  active  cracks  will  be 
dominant. 

An  approach  based  upon  rapid,  approximate  methods  for  determining  stress  intensity 
provides  good  indications  of  behaviour  provided  near-surface  stress  gradients  are  not 
excessive.  In  cases  where  stress  gradients  are  high  it  is  shown  that  it  is  necessary  to 
employ  numerical  techniques  in  calculating  stress  intensity. 

Two  specific  examples  are  presented,  the  first  at  very  small  scale  (heat-check  cracking, 
typical  crack  spacing  1mm)  and  the  second  very  large  scale  (permafrost  cracking, 
typical  crack  spacing  20m). 

The  predicted  ratios  for  the  proportion  of  cracks  shed  and  for  crack-spacing/cracktlepth 
are  in  agreement  with  experimental  evidence  for  gun  tubes,  concrete  and  permafrost. 
The  ratios  also  appear  to  match  experimental  observations  of  'island  delamination'  in 
ceramic  coatings  and  paint  films. 
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Figure  1 :  Surface  Topographies,  (a)  Dried-out  Mud  Flats  in  Death 
Valley,  typical  plate  size  250mm  (b)  Craze-Cracks  (Heat-Checking)  on 
Bore  of  a  Gun  Tube,  typical  plate  size  0.6mm  [3]. 
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(a) 


(b) 


(c) 


Figure  3 : 


(a)  Array  of  equaHength  edge  cracks 

(b)  Situation  after  first  shedding  of  2/3  of  cracks 

(c)  Situation  after  second  shedding  of  2/3  of  cracks 
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Figure  5  :  Variation  of  Stress  Intensity  with  h/a. 

Lower  curve,  h  fixed,  a  varying;  upper  curve,  a  fixed,  h  varying 
Filled  Squares  -  Situation  After  Further  Crack  Deepening. 
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Figure  8  :  Stress  Distribution  and  Associated  Approximate  Stress  Intensity  Factors  for 
Typical  Arctic  Permafrost  Scenario  in  Early  February. 
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Figure  9  :  Stress  Distribution  and  Associated  Accurate  Stress  Intensity  Factors  for 
Typical  Arctic  Permafrost  Scenario  in  Early  February.  (MMC  Solutions) 


Figure  10  :  Micrograph  of  Edge-Cracking  in  a  Gun  Tube. 
Depth  of  large  cracks  0.5mm  (after  [3]) 
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